Has Nucleonic Matter at the Critical Point Been Produced in Recent 
Multifragmentation Experiments? 
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For nuclei in five different mass regions the critical temperatures, at which the surface tension 
vanishes, are derived from information on the liquid branch of the coexistence curve. These critical 
temperatures increase with increasing mass and, for higher masses are well above recently reported 
critical temperatures obtained from Fisher Droplet model and percolation model analyses. However, 
for the lowest mass region nuclei, with 30 < A < 60, this analysis indicates that nucleonic matter 
has been produced essentially at the critical point, characterized by T c and p c . 

PACS numbers: 24.10.i,25.70.Gh 



The experimental and theoretical study of possible 
critical behavior in a nuclear fluid has occupied a consid- 
erable place in many contemporary studies (l], ||, j|, |], [|, 
[| 0, ^]. Given the complexity of nuclear collisions there 
is no a priori reason why collisional excitation of a nu- 
cleus, with its associated dynamic evolution should lead 
to a trajectory through the critical point of the nucleonic 
matter phase diagram and ascertaining whether the ac- 
tual critical point has been reached can be very difficult. 
Nevertheless, several papers have presented evidence for 
the observation of critical behavior in data from reactions 
of 1 GeV/nucleon 197 Au, 139 La and 86 Kr with 12 C, taken 
by the EOS group[[)| [II]] and data from reactions of high 
energy pionsp"l"|] and protons p~2| with 197 Au taken by the 
ISiS group. A Fisher Droplet model analysis (T3, Q of 
multifragmentation data for 8 GeV/c tt + + ™Au re- 
sulted in a critical temperature of 6.7 ± 0.2 MeV and 
p c « 0.3/Qo ji)- In a percolation analysis of 10.2 GeV/c 



197 



Au data critical temperature of 8.3 ± 0.2 MeV 
was determined pj. In a very recent paper 0, results 
from a Fisher Droplet model analysis of the nuclear "va- 
por" have been combined with a very interesting analysis 
based upon the Principle of Corresponding States, first 
elucidated for behavior of various atomic and molecu- 
lar gases |jl| In that work critical densities, p c — 
(0.39 ± 0.01)po and critical temperatures of 7.6 ± 0.2, 
7.8 ± 0.2, and 8.1 ± 0.2 MeV have been reported, respec- 
tively, for products of the reactions of 12 C with 197 Au, 
139 La and 86 Kr nuclei@. 

In this letter we apply the corresponding states analy- 
sis to information derived from the liquid branch of the 
coexistence curve. We conclude that the critical temper- 
atures of nuclei increase significantly with increasing nu- 
clear mass. The values we derive are generally well above 
the reported critical temperatures, except in the lowest 
mass range analyzed. In that mass range, 30 < A < 60, 
it appears that nuclei have been produced at densities 
and temperatures which are at, or very close to, those of 
the critical point. 

We have recently carried out analyses of existing 
caloric curve data obtained in a variety of nuclear re- 



action studies |L7||. These analyses provided evidence for 
mass dependent limiting temperatures and excitation en- 
ergies for nuclei. These limiting temperatures, and exci- 
tation energies are presented in Table 1. The temper- 
ature values were found to be in good agreement with 
predicted Coulomb instability temperatures |is|, [l9| , |2C|l 
calculated with theoretical models employing modern mi- 
croscopic nucleon nucleon interactions [|TJ |2^| . This gen- 
eral agreement has been exploited to derive the crit- 
ical temperature of symmetric nuclear matter, T c = 
16.6 ± 0.86 MeV as well as other information on the nu- 
clear equation of state p3|. 

The critical temperatures and excitation energies de- 
rived in references [ |, [| (?J are also listed in Table 1 to- 
gether with the masses of the disassembling nuclei actu- 
ally sampled at those critical excitation energies || [ll], 
O , [ijj . They are seen to be reasonably close to the lim- 
iting temperatures derived from the caloric curve mea- 
surements in the same mass range j^. Limiting temper- 
atures resulting from Coulomb instabilities would nor- 
mally be expected to fall below the critical temperatures 
of the corresponding nuclei as can be easily demonstrated 
in theoretical calculations by turning off the Coulomb 



TABLE I: Derived Temperatures and Excitation Energies 



Limiting and Critical Parameters (This work) 


Mass T Hm [17 


1 rj^n ucLeonic „ i t^* 
1 1 critical P/P0 & 


45 9.00 ± 1.17 8.52 ± 0.11 0.41 ± 0.004 7.71 ± 1.0 
80 6.97 ± 0.60 8.67 ± 0.92 0.42 ± 0.016 4.26 ± 0.5 
120 6.28 ± 0.50 10.1 ± 0.87 0.42 ± 0.009 3.58 ± 1.0 
160 6.63 ± 0.71 12.4 ± 0.99 0.41 ± 0.005 3.99 ± 0.5 
210 5.84 ± 0.41 12.8 ± 2.99 0.42 ± 0.015 2.53 ± 0.3 


Critical Parameters (Earlier work) 


Mass T C rit 


Reference Pc/po E* rit 


65 8.1 ±0.2 
106 7.8 ± 0.2 
140 7.6 ±0.2 
160 6.7 ±0.2 
168 8.3 ±0.2 


0.39 ±0.01 4.6 ±0.2 
0.39 ±0.01 4.9 ±0.2 
0.39 ±0.01 5.1 ±0.2 
0.3 3.8 ±0.3 
3.45 
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interaction|l8|. Thus the closeness of the temperatures 
at which apparent critical behavior is observed to the lim- 
iting temperatures seen in the caloric curves is somewhat 
surprising. In assessing the significance of the differences 
seen it should be noted that the excitation energies and 
the temperatures derived from the caloric curve analyses 
are those corresponding to the onset of the entry into 
a caloric curve plateau |l7j while the excitation energies 
extracted from the percolation and Fisher droplet model 
analyses are those corresponding to the points where ap- 
parent critical behavior is observed(5j [7]]. In those 
works the temperatures are then derived from the ex- 
citation energies assuming a Fermi gas behavior and as- 
signing a level density parameter. As a result these latter 
temperatures are not necessarily the same as the experi- 
mental values reported at those excitation energies. This 
by itself may be responsible for some of the tempera- 
ture differences observed in Table 1 and suggests that 
a comparison of the excitation energies is perhaps more 
meaningful. As seen in the table the critical excitation 
energies tend to be slightly higher than the excitation 
energies at the entrance into the plateau. 

In reference |24|] it was pointed out that an analysis 
of the caloric curves, carried out assuming a nondissi- 
pative uniform Fermi gas model, indicates a rapidly in- 
creasing expansion of the nuclei with increasing excita- 
tion energy above the excitation energy where the limit- 
ing temperatures are reached. In that work an iterative 
technique was employed to derive self consistent values of 
the thermal excitation energy, eth, and the relative den- 
sity, p eq /po- The densities derived in the early part of the 
plateau region of the caloric curve region were found to 
be w 0.60 — 0.75pof24|- Such densities are generally well 
above the values of 0.3 poto 0.39po derived at similar 
excitation energies from the Fisher Droplet and Corre- 
sponding States analyses || fjj. Only as the excitation 
energy is further increased are such densities reached in 
the caloric curve analysis. (See Table 1.) Thus it ap- 
pears that the results from the two different methods are 
in contradiction. 

To further explore this contradiction we have at- 
tempted to carry out a corresponding states analysis re- 
quiring that the results be consistent with our knowledge 
of temperature and density information for the liquid 
branch of the coexistence curve. The temperatures are 
readily available from the caloric curve but, for excitation 
energies below the plateau region of the caloric curve, the 
technique we have employed to derive densities is not 
applicable as shell, surface and effective mass effects re- 
duce the inverse level density parameter well below its 
Fermi gas value |^5|, |2(|. However, as pointed out in 
reference |24|] densities estimated from a model which as- 
sumes a trapezoidal density distribution and phenomeno- 
logical dependencies on temperature derived from previ- 
ous works are in good agreement with results of micro- 
scopic model calculations. Thus in the following, we es- 



timate the average densities at excitation energies below 
the onset of the plateau region using this model. Wc 
then take these average densities, together with the ex- 
perimental temperatures determined at those excitation 
energies as defining the liquid branch of the liquid vapor 
coexistence curve. 

As pointed out in references [|, 0, |7], the Principle 
of Corresponding States as developed by Guggenheim 
for inert gases results in a universal empirical scaling of 
Pi,v/Pc vs. T/T c of the form. 

Pl.v/Pc = 1 + 6i(l - T/T c ) ±b 2 (l- T/T c f (1) 

Where bi and b2= (1+ bi) are empirical constants, 
(3 is a constant equal to 0.33, the + sign between the 
second and third terms is valid for the liquid(l) branch 
and the - sign is valid for the vapor(v) branch of the 
coexistence curve. In the original scaling proposed by 
Guggenheim for macroscopic systems biis ss 0.75. El- 
liott et al. have pointed out that the constant appro- 
priate to a microscopic system is less as is demonstrated 
by with Ising model calculations for different sized cubic 
lattices 0. The analysis of Elliott et al. indicates a value 
of bi near 0.3 for the nuclei studied]?]]. In the following 
we have derived values of T c with fixed values of bi = 0.3 
and b2 = 1.3 and (3 — 0.333. We have also evaluated the 
uncertainties on the derived values of T c and p c which 
results from varying these parameters slightly. 

To illustrate the procedure which we have followed, we 
present in Figure la the results obtained by Elliott et 
al. from the analysis of the vapor branch obtained in 
1 GeV/nucleon heavy ion reactions H. This coexistence 
curve is constructed from the vapor branch assuming the 
critical temperature, T c , is equal to 7.6 MeV0. The va- 
por data (which are represented by triangles) shown in 
the figure are representative points taken from reference 
. These points fall exactly on the coexistence curve as 
they are the points from which the curve was derived. 
Also shown in the figure are the data points derived for 
the liquid branch in the in the 140 < A < 180 mass re- 
gion, using the approach described above (represented by 
circles) and the points in the caloric curve plateau, which 
are derived using the expanding Fermi gas hypothesis 
(represented by squares) (24| . It is immediately apparent 
that the liquid branch points do not fall on the coexis- 
tence curve. To fit such points a significant change in 
T c , p c , or both, is required. In Figure lb we present 
the result of fitting our liquid branch points of the coex- 
istence curve to equation 1. This fit leads to a critical 
temperature of 12.4 ± 0.99 MeV and a critical density of 
(0.41 ± .005)po- Thus the liquid branch implies a much 
higher critical temperature and the overall picture which 
emerges from such a fit indicates that the limit of the 
Coulomb instability prevents the system from reaching 
the critical point. Disassembly in the plateau region of 
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FIG. 1: Liquid gas coexistence curves. The solid line 
in each part of the figure depicts the liquid-gas coexistence 
curve determined in reference |7[ Symbols indicate data points 
from the vapor (triangles), liquid (circles) and mixed phase 
( squares) regions. The vapor points are representative data 
points from reference M. Liquid and mixed phase points are 
for nuclei with 140 < A < 180. See text. In part (a) the data 
are normalized to T c = 7.6 MeV and p c — 0.39po- In part 
(b) the critical temperature and critical density have been ad- 
justed to fit the liquid branch of the curve. There, T c = 12.4 
MeV and p c = 0.41po. 



the caloric curve would then appear to be occurring in 
the coexistence region well below the critical point. 

The results of extending this same method of analysis 
to data from the other mass regions are shown in Fig- 
ure 2(a)-(d). The temperatures and densities extracted 
from the fits to the liquid branches are presented in Ta- 
ble 1. While the critical densities remain near 0.4po, the 
critical temperatures increase with increasing mass. For 
all but the lightest systems (30 < A < 60) these critical 
temperatures are increasingly above those reported from 
the vapor branch analysis and well above the limiting 
temperatures derived from the caloric curves. 

However, for the lightest mass systems this analysis 
leads to a critical temperature of 8.5 ± 0.11 MeV which 
agrees quite reasonably with the 8.1 ± 0.2 MeV critical 
temperature extracted by Elliott et al. for their lightest 




P/Pc P/Pc 
FIG. 2: Liquid gas coexistence curves for nuclei in four ad- 
ditional mass regions. The solid line in each part of the figure 
depicts the liquidgas coexistence curve determined in refer- 
ence (?]. Symbols indicate data points from the vapor (trian- 
gles), liquid (circles) and mixed phase (squares) regions. See 
text. In each part of the figure the critical temperature and 
critical density have been adjusted to fit the liquid branch of 
the curve. See Table l.The Mass regions represented are (a) 
30 < A < 60, (b) 60 < A < 100, (c) 100 < A < 140 and (d) 
180 < A < 240. 



target, 86 KrQ . Here it should be recalled that the masses 
of the nuclei studied in that technique range from 80 to 
30 as the excitation energy deposition increases || [l0| |. 
The derived T c for 30 < A < 60 is also consistent with 
the limiting temperature of 9.0 ± 1.2 MeV obtained from 
the caloric curve analysis in this mass region ]i~7|]. Thus 
Figure 2(a) provides a strong indication that the criti- 
cal point may well have been reached in these light nu- 
clei where the Coulomb effects are less important. For 
the heavier systems the present results suggest that the 
apparent critical behavior observed actually occurs far 
from the critical point. This is consistent with the onset 
of important Coulomb instabilities |, ||, ||, ||. 

Recent theoretical investigations have indicated that fi- 
nite size effects or growth of fluctuations in the spinodal 
region may mimic critical behavior p7|, B8|. 

The results of this analysis arc summarized in Figure 3. 
There the critical temperatures derived from fits to the 
liquid branch are compared with the limiting Coulomb 
instability temperatures. Also indicated by a horizontal 
bar in the figure is the critical temperature for symmet- 
ric nuclear matter, previously derived from the limiting 
temperature data. While the limiting temperatures drop 
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FIG. 3: Mass dependence of limiting and critical temper- 
atures. The measured limiting temperatures (circles) and 
derived critical temperatures(squares) for nuclei are plotted 
against average mass for the different mass windows analyzed. 
The dot-dashed lines through the squares (circles) is a linear 
fit (2nd order polynomial fit) meant to guide the eye. Also 
indicated are the symmetric nuclear matter critical tempera- 
ture derived from these limiting temperatures (solid horizon- 
tal line)§| and the uncertainty associated with that value 
(dotted horizontal lines). 



with increasing mass, reflecting the increased importance 
of Coulomb effects, the critical temperatures rise toward 
the infinite nuclear matter value. The higher values for 
heavier nuclei are in accord with the constraint that T c 
should be greater than 10 MeV as suggested by Kar- 
naukhov based upon an analysis of fission probabilities 
taking into account the temperature dependence of the 
surface energy and its effect on the fission barrier |29| 
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